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The connection between the time-dependent physical spectrum of light and the phase space overlap 
of Wigner functions is investigated for optical pulses. Time and frequency properties of optical pulses 
with chirp are analyzed using the phase space Wigner and Ambiguity distribution functions. It is 
shown that optical pulses can exhibit interesting phenomena, very much reminiscent of quantum 
mechanical interference, quantum entanglement of wave packets, and quantum sub-Planck structures 
of the time and frequency phase space. 

I. INTRODUCTION 

In a paper published in 1977, Professor Eberly and one of the present authors (KW), have introduced the time- 
dependent physical spectrum of light described by non-stationary random electric fields The definition of such 
an operational spectrum required, as an essential step,a frequency tunable filter, that allows the resolution of the 
frequency components. If one restricts the result of |l| to deterministic fields only, the time-dependent physical 
spectrum takes the following form 



G(t,w) 



J dsJ{t-s-T)E{s)e %UJS 



(1) 



where E{t) is the positive frequency part of the detected electric field, and the spectral properties of the filter are 
represented by the response function J(t\ T)e lult . This function depends on the setting frequency ui, and the response 
amplitude is characterized by a bandwidth T of the filter. If one rewrites the time-dependent spectrum formula Q 
in the following form 



G(t,r) = 



(2) 



and look at as a function of time and bandwidth only, we recognize in this expression the modulus square of the 
wavelet transformation of the signal field, where t is the time shift and T is the scale 0- 

Several years after the introduction of the time-dependent spectrum, similar expressions to Q have been derived 
and applied to the reconstruction of the amplitude and phase of optical pulses. If one selects J{t) = \E{t)\ 2 or 
J{t) = E(t), the recorded time-dependent spectrum corresponds to the frequency resolved optical gating (FROG) 
0, or to the second harmonic frequency resolved optical gating (SHG FROG) Methods based on FROG have 
become powerful tools in the investigations of femtosecond pulses. 

The remarkably simple expression Q , hides a very interesting phase space structure of the operational spectrum 
in terms of a time and frequency distribution overlap. It has been shown p|, that the time-dependent spectrum can 
be visualized as a time and frequency convolution of two Wigner distribution functions 

G(t,u) =2n J d£ J dui' 'Wj(t' \u')W E {t' +t,w' + u) . (3) 

In this expression the first Wigner function corresponds to the time and frequency distribution of the filter, while 
the second Wigner function describes the phase space properties of the measured electric field. It is perhaps worth 
mentioning, that a quantum version of the expression J2J, can be applied in quantum mechanics, to describe joint 
operational measurements of position and momentum |y| . 

In order to honor Professor Eberly's contributions to the development of the time-dependent spectrum and Quantum 
Optics, this paper will investigate time and frequency properties of optical pulses using the concept of the phase space 
Wigner distribution. We will show that optical pulses can exhibit interesting phase space structures very much 
reminiscent of quantum mechanical interference 7], quantum entanglement of wave packets |8(, and quantum sub- 
Planck structures [!j ■ 
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This paper is organized in the following way. Section II is devoted to the definition and elementary properties of 
the time and frequency Wigner and Ambiguity distribution functions. Section III contains a detailed description of 
chirped pulses using a time and frequency phase space. We show that time and frequency correlations of chirped 
pulses have a formal analogy to entanglement of wave packets. The strength of these correlations is investigated using 
the Schmidt decomposition. In Section IV we analyze the phase space properties of a linear superposition of two 
chirped pulses. In Section V we discuss the connection of sub-Fourier phase space structures with pulse overlaps. 
Finally some concluding results are presented. 



II. TIME AND FREQUENCY PHASE SPACE 
A. The Wigner Function 

We shall investigate time and frequency properties of optical pulses using a phase space distribution function. Such 
a function has been originally introduced by Wigner in 1932 and applied to quantum mechanics jlCj . In the area of 
signal processing the same distribution function has been used by Ville in 1948 |ll|. The time and frequency Wigner 
distribution function corresponding to the field envelop E(t) is defined as 

W, M = J #>-E(t-$. (4) 

It is well known that the Wigner function can be used as a time and frequency distribution of the pulse, but cannot 
be guaranteed to be positive for all fields. There is an extensive literature devoted to the properties and applications 
of the Wigner function in quantum mechanics [T^ and classical optics |l3j | . Below we present only the most relevant 
properties of the Wigner distribution needed for the purpose of this paper. 

The frequency integration of the Wigner function yields the temporal instantaneous intensity 



dujW E {t,u) = I(t) = E*{t)E{t) . (5) 
The corresponding time integration of this distribution leads to the power spectrum of the optical pulse 

/ <toW E (t,u) = P(u) = ^E*(u)E(w). (6) 

J 27T 

In this formula the expression E(u>) = J dte~ E(t) is the Fourier transform of the pulse. From these definitions we 
see that that the Wigner function is normalized to the total power/energy of the pulse 

J dt J duW E {t,uj) = J dtl(t) = J duP(uj), (7) 

where the last equality follows from the Parseval theorem for the Fourier transforms. An important result that we 
shall use in the following sections is the overlap relation for two Wigner functions 
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\{Ei\E 2 



dtE*(t)E 2 (t) 



2tt / dt / dwW El (t,u)W E2 (t,u). (8) 



This formula indicates that the case a zero overlap of two pulses is impossible to achieve with positive Wigner functions. 

Using the Wigner function as a weighting distribution one can characterize the properties of the optical pulse in 
the form of the following statistical moments of time and frequency 

[ ' ~ JdtJdtoW E (t,Lu) ' [ ' 



B. The Ambiguity function 

A different way of looking at the time-frequency correlations © is to use the Ambiguity function, which is a 
two-dimensional Fourier transform of the Wigner function 



A E (T,Q) = j dt J duje int+tujT W E (t,u) 



= I <HE ! (*-|) e^Eft+j) (10! 
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From this definition it follows that the Ambiguity function can be written as 

A E (T, O) = A E (0, o)(e* nt +- T ) , (11) 

where A E (0, 0) = J dt J dio W(t, ui) is a normalization constant. The formula (|11J) can be used as a moment generating 
function. The time and frequency statistical moments |j!|J) can be calculated from the Ambiguity function using the 
formula 



(12) 

T=n=o 



C. ABCD optics of optical pulses 

It is known from classical optics that for linear optical devices one can use the ABCD transformation of geometrical 
ray displacement and slope [l4j | . We will use this approach to describe arbitrary linear transformations of time and 
frequency given by 

50 -(s !)(:)■ ™ 

This transformation is canonical if it preserves the normalization of the Wigner function 

J dt J dwW(At + Bu,Ct + Duj) = J dt J dioW(t,Lu). (14) 



detl (I.r,i 



This condition is satisfied if 

Ac+ ( 

C* D 

The ABCD transformation of the Wigner function generates the following transformation of the Ambiguity function 

W{At + Bui, Ct + Duj) =*> A E {At - Bui, Dui ~ Ct) . (16) 

III. TIME AND FREQUENCY DESCRIPTION OF CHIRPED PULSES 
A. General properties of chirped pulses 

Let us write the electric field in the form of a real envelop and phase: 

E{t) = £{t)e~ lv ^ . (17) 

In order to calculating the Wigner function for such a pulse we use a linear approximation for the phase: ip(t + s/2) — 
tp(t — s/2) ~ ips, and as a result we obtain 

W E (t,u) = W e (t,u + (p), (18) 

where Ws {t,u>) is the Wigner function of the real envelop £ (t) . The corresponding formula for the Ambiguity function 
is 

A E (t,u) = Ae(t,u-(p). (19) 
The instantaneous pulse frequency ui{t) of a chirped pulse E(t) = £(t)e luJ ^\ can be defined as 

, , ( du> us W E (t, oj) , , . 
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and the corresponding square of instantaneous pulse frequency is 

where and (w 2 )f are the frequency two moments, calculated with respect to the unchirped amplitude £. From 
these relations we conclude that the dispersion of the instantaneous pulse frequency is 

(Au,(t)) 2 = (Acj)l + (Amf ■ (22) 



B. Gaussian optical pulses with linear chirp 

As an example of the general envelop l|17fl . we will consider a single pulse with a linear chirp and Gaussian envelop 
function of the form 

E{t) = cxp (-iuit - JL(l + ia)j . (23) 

We have assumed that our pulse is long enough, so one can perform the standard decomposition of the electric held 
into the slow amplitude (|17|l . and the harmonic carrier with frequency u>i. The intensity and the power spectrum of 
this pulse are 

/ t 2 \ , . 2a 2 ( 2cj 2 a 2 \ 

/(t) = «p(-^j, PM = - 7 _exp(- IT ^j. (24) 

In these formulas the full duration of the pulse is defined as a full width at half maximum of the intensity (FWHM) : 
t p = 2\/2 In 2 <7, the linear chirp is characterized by a real parameter a, and the electric held amplitude has been 
conveniently selected to be one in arbitrary units. In all numerical applications in this paper, we select a 2 = 1/2 
leading to a pulse duration 2-s/ln 2. In order to keep our formulas simple we have shifted the frequency u> in such a 
way that it incorporates the constant carrier frequency. 

The chirp on the pulse (IL'MI) corresponds to a linear chirp tp(t) = i&st, leading to the instantaneous pulse frequency 
w(t) — —^it- For the Gaussian pulse the formula Q22|l becomes 

( A -) 2 = i + 2^ 2 >- ( 25 ) 
We see that the linear chirp is equivalent to a transformation given by the following ABCD matrix 

(26) 

This matrix has the typical form for a ray transformation due to a thin lens |14| . 

Let us investigate the time and frequency properties of the chirped pulse. Simple calculation shows that the Wigner 
function of this pulse is 

W(t,w) = W (t,u+^t) 

= -^=expf — ^— (1 + a 2 ) - 2uj 2 a 2 - 2aut] , (27) 
V2^ V 2cr / 

where Wo(t, w) — -j== cxp (^— — 2uj 2 a 2 ^j is the Wigner function of a Gaussian pulse with no chirp. In Figures Q 

and J5J we have depicted the Wigner function of a Gaussian pulse with no chirp and chirp a = 3. 

For the Ambiguity function a simple calculation for the chirped pulse gives an expression similar to the formula 

m 

A B {T,n) = A E (0,0)A o (T,n- ^T) 

= A E (0,0)ex P f-^(l + a 2 )-^ + ^Tn\ . (28) 





2-1 


f 1 








a 

V 2a 1 





5 




This function can be written as a general Gaussian function of two variables 

1 



A E (T,n) = A £ (0,0)exp 



Cv T 



(29) 



where v = (T, Q) is a vector and C is a 2 x 2 covariance matrix of the time and frequency variables. As a result 



(t 2 ) <M 

(iw) (u, 2 ) 



a 1+q 2 
" 2 4cr 2 



Note that 



detC 



From this covariance function we obtain that time and frequency dispersions are 



At = a, Aw = — \J\ + a 2 , 
2a 

and that the Fourier uncertainty relation between frequency dispersion and time dispersion is 

. . Vl + a 2 1 

At Auj = — — > - , 

2 ~ 2 ' 



(30) 



(31) 



(32) 



(33) 



where the lower bound corresponds to Gaussian pulses with no chirp. From these relations we see that the chirp 
enlarge the spectral width of the pulse. 
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From the covariance matrix (|30|l we see that the chirped pulse leads to the following time-frequency correlation 

{tu) = ~. (34) 

In the next Section we show that the Ambiguity function is particularly useful to quantify the "strength" of this 
time-frequency correlation. 



C. Schmidt's decomposition of chirped pulses 

We will use the Schmidt decomposition in order to quantify the correlation properties of the chirped pulses. This 
decomposition has been successfully used to quantify entanglement of quantum mechanical systems described by a 
correlated two-party wave functions. As an example we note that recently Professor Eberly has investigated high 
transverse entanglement in optical parametric down conversion using a Schmidt decomposition of the biphoton wave 
function 

In order to apply the Schmidt decomposition we replace (T, ft) by two dimensionless variables (X, Y) such that the 
Ambiguity function takes the following form 

A(X, Y) = A(0, 0) exp(-A 2 — Y 2 — 2cXY) , (35) 

where 

c=^t^. (36) 

The Schmidt decomposition of the function (fTUfl is 



A(X,Y) = A(0,Q)J2VKMX) 1 Pn(Y) (37) 

n 

where ip n (X) and ip n (Y) are the Schmidt modes defined as eigenstates of the reduced density operators constructed 
from a two-party wave function given by the Ambiguity function (|35|> . The Schmidt eigenvalues p n serve as a degree 
of entanglement or correlation between the two-part system. In our case the Schmidt eigenvalues quantify the degree 
of correlation between time and frequency. The measure of correlation can be quantified by entropy of entanglement 

E = l-^pl =*> Q<E<1. (38) 

n 

For a untangled pulse we have E = because there is only one Schmidt eigenvalue, and the formula l|37() factorizes: 
A(X, Y) = A(X)A(Y). 

The Schmidt decomposition can be performed for a Gaussian ambiguity function characterized by an arbitrary 
covariance matrix (|3(J|I . In this case the entropy of entanglement is 



E=i -f~M^' (39) 

Clearly we have E = for uncorrelated in time and frequency optical pulses. 

For the Gaussian optical pulse (|23|l with chirp, we obtain that the Schmidt eigenvalues have the formal form of the 
Bose-Einstein distribution 



Pn 



(l + n)™+ 1 ' ' "' 2 



with n = -(y/l + a 2 - 1) . (40) 



The Schmidt eigenstates corresponds to a thermal distribution with a mean number of photons given by n. Using 
this expression we calculate the entropy of entanglement (|39|l 

E = 1 — - 1 . (41) 



In Figure © , we have depicted the entropy of entanglement as a function of the chirp a. 
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We conclude this section with a remark that pulses with linear chirp are analogous to quantum mechanical wave 
functions of the form ip(x) = e lax . In the framework of quantum mechanics such wave functions belong to a wide 
class of the so-called contractive states [l^, and have been used to show the narrowing of the position uncertainty 
relation , and applied by Professor Eberly to phase entanglement ^3] . The contractive nature of the chirped pulse 
can be easily exhibited using an lens type ABCD transformation of the chirped pulse, equivalent to the following 
linear transformation in time and frequency plane 

u t = w + ■ (42) 
For a chirped pulse the dispersion of this observable is 

(A Wt ) 2 = (Ac) 2 + -^(Atf + \{u>t) = 2 -±^-«. (43) 

For positive a, this formula exhibits a contraction of the uncertainty of u>t similar to the narrowing of a freely evolving 
quantum mechanical wave packets. The narrowing is entirely due to the time-frequency correlation of the chirped 
pulse JjUJl. 



D. Time dependent spectrum for chirped pulses 



In order to calculate the time-dependent spectrum for the chirped pulse, we will use for the filter response 



4r 5 



J{t)=0(t)T\ -exp 



(44) 



This function corresponds to a linear, causal and time translation invariant response of the filter. In the calculations 
we will assume that the observation time is larger compared with the pulse duration: t > T4. Simple calculation shows 
that the spectrum is Gaussian and has the form 



G(t, uo) — G exp 



1 



w G w T 



(45) 



where v = (t, ui) is a vector and the matrix G can be expressed in terms of the covariance matrices of the pulse and 
the filter 



G 



C + Cj 



1 1 l+a^ 



The normalization constant Go of the time-dependent spectrum is such that 

dtdbj 



2tt 



G(t,tu) 



dt I(t) 



(46) 



(47) 
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IV. LINEAR SUPERPOSITION OF CHIRPED PULSES 



Linear superposition principle play a fundamental role in classical and quantum interference phenomena. In order 
to illustrate time and frequency interference we shall investigate a linear superposition of two optical pulses 

E sup = E(t - t ) + E(t + t Q ) , (48) 

where 2to is a temporal separation between the pulses. This linear superposition of two electric fields of optical pulses 
exhibit classical interference very similar to the interference of quantum coherent states [Isl Il9l | . 
As an example of such a superposition we take two Gaussian chirped pulses 

S sup (*)=e X p('-^^(l + m))+exp('-^i^(l + m)) . (49) 

The intensity of this superposition is: 

The corresponding power spectrum is 



4(T 2 / 2uo 2 a 2 \ 2/ , 

P(oj) = - 7 === cxp -— - cos 2 Mo • 51 
V 1 + a 2 V 1 + a J 

The time-frequency Wigner function of the linear superposition Q49J1 is 

W{t,uj) = W E {t + t ,uj) + W E {t-t Q7 uj) + 2W E {t,uj)cos(2t uj). (52) 

In this formula the Wigner function: We^, u>) is given by H27|) . From the Wigner function of the superposition it is 
possible to calculate time and frequency moments. Simple calculation leads to 

tg(l-a 2 exp(-%^)) 

(* ) - a + — : . ,j^2,, — , ( 53 ) 



1 + exp(— — 



2a- 



and 



1 + a 2 ( t 2 2 exp( 



*5(i+" )' 



4a V l + exp(-^^i)/ 

From this formula we see that the spectrum of the linear superposition is reduced {squeezed) below the single pulse 
width. This effect is entirely due to the fact that we are dealing with a linear superposition of pulses. In Figures Q 
and |J5J we have depicted the Wigner function for the linear superposition with to = 4 with no chirp and chirp a = 3. 
The squeezing effect corresponding to the nonzero chirp is clearly seen. 

Note that in reference |2(| , the Wigner function of a coherent two- pulse sequence with linear frequency chirp has 
been reconstructed experimentally using quantum tomography. 

The remarkable feature of the Wigner function l|52ll is the fact that it contains structures in a phase space region 
below the Fourier uncertainty relation. In Figure © we have depicted the Wigner function in a space region with 
Auj < 1/ V2 and At < 1/ y/2. In the framework of quantum mechanics it has been recognized that small structures 
on the sub-Planck scale do show up in quantum linear superpositions 0- It is clear that for linear superpositions of 
chirped pulses such sub-Fourier structures emerge as well. We will see in the next Section, that due to such small 
structures it will be possible to have pulses with zero overlap. 



V. TIME AND FREQUENCY OVERLAP 

In this Section we will investigate a FROG version (with J{t) = E* up (t)) of the time dependent spectrum Q 
ipplied to the linear superposition of pulses given by the formula 1491) 



G(0,A) = 



J ds E; up (s)E sup (s)e % 



(55) 
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FIG. 4: Plot of the Wigner function for two Gaussian pulses with no chirp a = and to = 4. 




FIG. 5: Plot of the Wigner function for two Gaussian pulses with chirp a = 3 and to = 4. 



Such an overlap can be easily obtained if one of the two pulses has its carrier frequency detuned by A. In this case 
the FROG signal is a phase space overlap JSJ) of the form 

G FROG (0, A) = 2irjdtj diu W Eenv (t, lo)W Es ^ (t, to + A) . (56) 

For the chirped pulses this overlap can be calculated, and as a result we obtain 

G FROG (0, A) = Afe-^~ I cos(Aio) + e _ ^ (1+ ° } cosh(Ai a) J , (57) 

where A/" is a normalization constant. In Figure Q we have depicted l|57|l as a function of A, for to = 5 and a — 5. 

The remarkable feature of this curve is that GV.rog(0, A) can vanish for values of A smaller then the frequency 
uncertainty in frequency. In terms of the time and frequency distribution function, the detuning A of the carrier 
frequency corresponds just to a shift of the Wigner function along the u> axis. In the case of the superposition of two 
Gaussian pulses, this shift affects the Gaussian peaks and the oscillating non-positive interference term. It is easy to 
notice that for an appropriate value of A this can result in a sign change of the interference term in comparison to the 
original function. Calculating the overlap (|56|) we need a product of shifted and unshifted Wigner functions. For a 
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FIG. 6: The interference terms of the Wigner function for two Gaussian pulses plotted at the region of the phase space below 
the Fourier uncertainty: a) fragment of the plot from Figure @, i.e. a = 0, to — 4; b) fragment of the plot from Figure JKJ, i.e. 
a — 3, to = 4; c), d) similar terms obtained for larger time separation distance to — 9 with chirps a = and a = 5, respectively. 



G FROG (0,A) 




1 2 3 4 5 A 
FIG. 7: The overlap as a function of A for a linear superposition with to = 5 and a — 5 



selected value of the shift A, the product of the two interference terms of their Wigner functions will become negative 
(or zero). In Figure JHJ we have depicted such a product with to = 4 and shift A = ir/2t . We see that the considerably 
large negative contribution can cancel the positive peaks corresponding to the overlap of the non-interference terms. 
Obviously, formula l|57f) can be easily obtained without using the Wigner function, but the phase-space representation 
(15611 is especially useful to show that the sub-Fourier structures are peculiar to interference phenomena. 

In quantum mechanical framework this feature is especially interesting as the zero-overlap means that the corre- 
sponding states are orthogonal, and so at a sub-Planck scale one can obtain sets of mutually orthogonal states 0. In 
the case of optical pulses we can obtain a similar result of pulses with zero-overlap using a shifts of frequency that are 
below the Fourier uncertainty relation. Certainly, the smaller is the value of this shift A that one wants to use, the 
larger to has to be taken in the calculations, meaning that the Fourier/Heisenberg uncertainty relation is satisfied. 
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FIG. 8: The product of two Wigner functions for linear superpositions of pulses with to — 4 and a — mutually shifted in 
carrier frequencies by A = 



VI. CONCLUSIONS 



In this paper we have analyzed the connection of the time and frequency time-dependent spectrum with phase space 
distributions based on the Wigner and the Ambiguity functions. We have exploited the similarities between optical 
pulses and quantum mechanical wave packets to exhibit interference, entanglement, and sub- Fourier structures of the 
time and frequency phase space. 

Certainly the Wigner function is not the only possible way to give a time and frequency description of optical 
pulses. Recently in reference |2l| we have introduced a class of new phase space representations based on the so called 
Kirkwood-Rihaczek distribution function : 



K(t,v) = 



-t) 



Eit) 



(58) 



Such phase space distributions provide time and frequency characteristics of optical pulses, but only phase space 
overlaps of these distributions have an operational meaning as a physical spectrum of light. In the case of the 
Kirkwood-Rihaczek function this overlap, much resembling the equation (JSJ, is given by the following formula 



dtE*(t)E 2 (t) 



= 2tt / dt / duK El (t i u))KZ a (t,u). 



Thus, for the superpositions discussed in Section V, we can write the following formula 

G FflOG (0,A) = 2ir J dt J duK Eeap (t,co)K* Emp (t,uj + A). 



(59) 



(60) 



We conclude this paper by noting that the definition of the time-dependent spectrum of light introduced 28 years ago, 
is still an attractive field of research producing interesting insights into the definition of the spectrum of nonstationary 
ensemble of optical pulses [2^ • 
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